Enhanced signal of astrophysical tau neutrinos propagating through Earth 
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Earth absorbs Ve. and of energies above about 100 TeV. As is well-known, although v T will 
also disappear through charged-current interactions, the v T flux will be regenerated by prompt tau 
decays. We show that this process also produces relatively large fluxes of secondary i/ e and £" M , 
greatly enhancing the detectability of the initial v T . This is particularly important because at these 
energies v T is a significant fraction of the expected astrophysical neutrino flux, and only a tiny 
portion of the atmospheric neutrino flux. 
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Introduction. — A variety of astrophysical and exotic 
sources are expected to produce large fluxes of ultra-high- 
energy particles, including photons, protons, and neutri- 
nos. Sources such as active galactic nuclei and gamma- 
ray bursts are at high redshifts, and the photons and pro- 
tons will be attenuated by scattering on cosmic radiation 
backgrounds. The protons may also be deflected by mag- 
netic fields. However, neutrinos from distant sources will 
neither be attenuated nor deflected, thus allowing true 
neutrino astronomy of the high-redshift universe. Fur- 
ther, it is now known that about one-third of the energy 
density of the universe is accounted for by particles out- 
side the Standard Model. Among the possibilities are 
variants of superheavy dark matter, which may produce 
ultra-high-energy particles from their annihilations or de- 
cays. Such signals would be most readily observed from 
local sources. The existence of cosmic rays with ener- 
gies up to and above 10 20 eV is a powerful argument for 
the existence of ultra-high-cncrgy astrophysical neutri- 
nos. With the upcoming km-scale detectors, a new era 
of neutrino astronomy will begin, and the energies, direc- 
tions, and flavors of astrophysical neutrinos will provide 
important clues to the most violent astrophysical objects 
and the nature of particle dark matter jl],|| . 

While these detectors can suppress downgoing atmo- 
spheric muons, atmospheric neutrinos are a much more 
challenging background. Like astrophysical neutrinos, 
these can pass through Earth and create upgoing muons 
(for example). They are produced in the decays of sec- 
ondaries produced by the collisions of cosmic-ray protons 
with Earth's atmosphere. The decays tt + /K + — > [i^v^ 
and — > e + h , e 9 fl and their charge conjugates produce 
fluxes in the flavor ratio v e : : v T ~ 1 : 2 : 0, where 
neutrinos and antineutrinos have been combined. These 
ratios vary with energy and angle when the production 
of other mesons and their decay lengths relative to the 
slant depth of the atmosphere are taken into account. 
In particular, at the highest neutrino energies, there is 
a tiny relative v T flux from charm decay The at- 
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mospheric neutrino flux falls roughly as E~ 3 , whereas 
general considerations for astrophysical neutrino fluxes 
predict shallower slopes, and hence that they will domi- 
nate above some energy. For diffuse fluxes, this crossover 
is thought to be at least 10 TeV (point and/or transient 
sources can be identified at lower energies) Jl],|2| . 

Neutrino oscillations have an important consequence 
for suppressing the atmospheric neutrino background. 
The Super-Kamiokande atmospheric neutrino results Q| 
strongly favor <-» v T oscillations with sin 2 28 ~ 1 and 
Am 2 ~ 3 x 10~ 3 eV 2 . At the high energies considered 
here, <-> v T oscillations will never occur for atmo- 
spheric neutrinos, and their flavor ratios will reflect the 
production mechanism only. In the astrophysical pro- 
duction scenarios, the flavor ratios are also thought to be 
1:2:0 (again combining neutrinos and antineutrinos), 
reflecting their production by ultra-high-energy protons. 
But for the astrophysical sources, <-> v T mixing will al- 
ways be complete because of the long path lengths, and 
so these ratios become 1:1:1. In some exotic mod- 
els [HJo), the initial v T flux is large; in any case, a large 
v r component in the astrophysical neutrino flux is guar- 
anteed. 

This simple fact is very important for the detection of 
astrophysical neutrinos above the atmospheric neutrino 
background. This is because v e and above about 100 
TeV will be absorbed by Earth, whereas the absorbed v T 
flux will be regenerated by prompt tau decays jjj. We 
point out a new effect: that the regeneration of the v T flux 
also creates relatively large secondary fluxes of v e and . 
These make the detection of astrophysical sources signif- 
icantly easier, and must be accounted for to properly 
deduce the source spectra and flavor composition. 

Neutrino propagation in Earth. — Above about 100 
TeV, all neutrino flavors have a high probability of in- 
teracting in Earth via neutrino-nucleon scattering, with 
a total cross section a oc E - 5 ||. Except for v e + 
er — > W~ at the Glashow resonance (6.3 PeV), neutrino- 
electron scattering is irrelevant. Charged-current (CC) 



1 



and neutral-current (NC) interactions on nucleons occur 
in the ratio 0.71 : 0.29, and in both cases, the surviving 
lepton carries about 75% of the initial neutrino energy ||. 
The cross sections for neutrinos and antineutrinos are 
approximately equal. The neutrino interaction length is 
shown in Fig. [IJ it is the same for all flavors since even 
tau mass threshold effects are negligible. 

In the more likely CC interaction, the critical ques- 
tion is the subsequent fate of the charged lepton. An 
initial v e will be removed from the beam and the result- 
ing electron very quickly brought to rest. With an initial 
v^, the produced muon will eventually decay. However, 
the laboratory-frame decay length for the muon is always 
much longer than its range, the distance over which its 
electromagnetic energy losses bring it to rest. Thus the 
neutrinos from its decay will be below 50 MeV, and of 
no further interest here. The decay length and range for 
muons are shown in Fig. [|. 

As first noted by Ritz and Seckel ||, the situation is 
very different for v Tl since at all but the highest energies, 
the tau decay length is less than its range, so that the 
tau will decay in flight without significant energy loss; 
see Fig. |l|. Tau decays occur by many branches, but all 
contain a v T . This regenerated v T carries a fraction about 
0.75 x 0.4 ~ 1/3 of the initial v T energy, where 0.4 is for 
v T from tau decay JT3JO]. Ritz and Seckel considered 
WIMP annihilation in the sun, and the propagation of 
the produced neutrinos. The utility of this regeneration 
effect for the detection of astrophysical neutrinos was not 
noted until a decade later, by Halzen and Saltzberg 0. 

The regeneration process will continue (along with the 
less frequent NC scatterings) until the v T energy has been 
moderated to the energy such that the neutrino interac- 
tion length is comparable to the remaining distance in 
Earth. At nadir angle ip = Q°, such that the neutrino 
will cross the full diameter of Earth, the transparency 
energy is about 40 TeV (at larger nadir angles, neutrinos 
at higher energies can pass through without scattering). 
The distribution of v T energies around the transparency 
energy is characteristically lognormal, with a one-sigma 
width of approximately one decade. Thus all of the initial 
v r will emerge at relatively high energies. 

Future km-scale astrophysical neutrino detectors are 
designed for the detection of v e (via an electromagnetic 
shower) or (via the track of a penetrating muon). 
The detection of v T is much more difficult. Learned and 
Pakvasa jl2] have shown that v T may be detected by 
a "double-bang" signature, if the tau lepton production 
and decay, each accompanied by a large shower, are well- 
separated but both occur in the detector. The double- 
bang events will be observable only in a narrow range 
around several PeV, and only at large nadir angles, so 
other techniques for v T detection are needed. 

This is precisely the point of the signal proposed by 
Halzen and Saltzberg Q. In their scheme, the v T inter- 
acts outside the detector and the tau created decays as 
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FIG. 1. Some important length scales, in units of Earth 
radii and km, versus the energy in GeV for the neutrino or 
charged lepton. Solid lines: the fi and r decay lengths, i.e., 
the distances over which a fraction e _1 of the the initial fluxes 
would survive (ignoring energy loss) . Dashed lines: the /x and 
t ranges, i.e., the distances in standard rock (22 g mol -1 ) at 
a density 8 g cm -3 over which the which they would be fully 
stopped by their electromagnetic interactions (ignoring de- 
cays). Dotted line: the neutrino interaction length (the same 
for all flavors), i.e., the distance at a density 8 g cm -3 over 
which a fraction e" 1 of the initial flux would survive. With 
the exception of the Glashow resonance for £> e at 6.3 PeV, this 
figure is the same for antiparticles. While the density varies 
significantly with nadir angle, 8 g cm" 3 is representative. 

t~ — > iiTVpVT, where the muon, which has a long range, 
is seen in the detector. These v T events can be sepa- 
rated from CC events (which also produce a muon) 
on a statistical basis, using the fact that the v T nadir 
angle distribution will be unchanged by passage through 
Earth, whereas that for will show exponential absorp- 
tion (as will v e ). Additionally, if the incoming neutrino 
spectrum is shallower than about E~ 2 , the pileup near 
the transparency energy may be visible [p]Jl0|. However, 
it should be noted that the above tau branching ratio 
of 18% reduces the detectability of the v T flux compared 
to a i/fj, flux at the same energy, both taken just below 
the detector. While the muon in the ^v-induced case 
will have an energy roughly 0.4 of that in the ^-induced 
case because of the tau decay |l(],|ll[], the muon range 
and hence detectability varies only slowly with energy 
(see FigQ). 

Secondary antineutrinos. — While expected to be a sig- 
nificant component of the astrophysical neutrino flux, 
and the only component which is not attenuated in 
Earth, v T are not easily detected. However, we point out 
a new effect that significantly improves their detectabil- 
ity. In each step of the regeneration process v r — > r — > v T 
described above, the tau lepton decay will always pro- 
duce a v T . Additionally, 18% of decays are r — > v-rtw^ 
and 18% are r — > v T ev e . In the laboratory frame, on av- 
erage v T and \x each carry a fraction 0.4 of the tau energy 
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and carries a fraction 0.2 JjjjjrjJ . Secondary and 
v e are also produced by initial D T (expected to have a 
flux comparable to the initial v T ). These secondary an- 
tineutrinos (p e , P M ) and neutrinos (v e , u^) have not been 
taken into account before, but they have a surprisingly 
large effect on the detectability of v T and v T . 

The number of regeneration steps for v r is approxi- 
mately N — \og(Ei / Et) / \ogi : where Ei is the initial 
energy and Et is the transparency energy. With more 
steps, there are more opportunities to create secondary 
and v e . However, those created farther from the de- 
tector have a greater chance of being absorbed. On each 
step the secondary antineutrinos will carry about 1 /6 of 
the initial v T energy, reflecting the energy lost in the CC 
scattering and the tau decay. 

It is generally assumed that all scatterings and decays 
are collinear. However, as the number of regeneration 
steps increases, the cumulative non-collinearity could be 
large enough to blur astrophysical point sources. If the 
initial neutrino has energy E v , simple considerations in- 
dicate that for v T — > r, the scattering angle is about 
1°/ \jE v j TeV, and for r — > v T , the decay angle is about 
0.2° l(E v I TeV). The decay angle can thus be neglected 
compared to the scattering angle. In order for the num- 
ber of scatterings to be nonzero, E v must be above about 
100 TeV, so that on the last step the angular deviation 
is about 0.1°. Taking previous regeneration steps into 
account gives a maximum deviation of about 0.3°. This 
is below the expected 1° reconstruction resolution of the 
proposed detectors 0] , and so astrophysical point sources 
should not be blurred by v T regeneration. 

We can estimate the secondary v e and number flux 
produced by a v r beam by making the following crude 
assumptions: (i) v r do not lose energy in their CC or 
NC interactions; (ii) v e or which interact are removed 
from the beam; and (iii) v e or from tau decay are 
produced with the same energy as the initial v T . Consider 
a beam of astrophysical neutrinos in the expected ratio 
v e : : v T = 1 : 1 : 1, each with a flux Fq, and all at an 
energy E . Define / as the fraction of v T interactions that 
produce a a v e or v^: f = 0.71 x 0.18 = 0.13, where the 
first factor is the CC/NC fraction, and the second factor 
is the relevant tau-decay branching ratio. The secondary 
i/fj, flux from tau decay obeys the equation 

where x is the distance through Earth at the chosen nadir 
angle. The first term represents the creation of from 
the v r flux, and the second represents their subsequent 
absorption. Note that the interaction length A is the 
same for both flavors. Thus, in this simple estimate, the 
fluxes after propagation through Earth are 

FujF = F V JF = e~*/ x ; F V JF = 1 

F„JF = F P JF = f (l - e-* /A ) ; F,JF Q = . 



The initial v e and fluxes are exponentially depleted. 
However, the flux associated with the initial v T is in- 
creased due to the secondary v e and 9^. For large x, cor- 
responding to more v T regeneration steps, the increased 
production of v e and is balanced by their increased 
absorption. For the expected initial antineutrinos (with 
the same fluxes and flavor ratios as for neutrinos), similar 
results obtain for the production of secondary neutrinos, 
and they must be added to the above. For small x, the 
initial v e , v^, D ei P M fluxes will survive, whereas for large 
x, they will be absorbed but replaced by the approxi- 
mately constant secondary fluxes created by v T and v T . 
Thus even for large energies and/or small nadir angles, 
the fluxes of v^, D e , never vanish, contrary to pre- 
vious predictions. 

We have calculated the secondary v e and produc- 
tion with a Monte Carlo code that simulates the pas- 
sage of high-energy neutrinos through Earth. We assume 
collinear propagation for all particles followed, as justified 
above. The code starts with an initial neutrino and ran- 
domly samples its interaction probability at each step in 
distance, using the CC and NC total cross sections from 
Gandhi, Quigg, Reno, and Sarcevic ||. In each interac- 
tion, we sample the outgoing charged lepton or neutrino 
energy around the average inelasticity values (y) given in 
Ref. B. At these energies, (y) is the same for CC or NC 
interactions, and the same for neutrinos or antineutrinos. 
Accordingly, for the distributions about (y) we use the 
da I dy distribution given in Ref. S for neutrino-nucleon 
CC scattering. 

After a NC interaction, the neutrino is followed un- 
til the next interaction. After a CC interaction, if a 
muon or tau is created, it is followed. For their energy 
loss, we use the results of Ref. |R|, and for their de- 
cays we use the laboratory-frame decay probability dis- 
tributions. Any electrons created are assumed stopped. 
For the laboratory-frame v r distribution from tau decays, 
summed over all branches and including the tau polar- 
ization, we use the distributions given in Ref. |ic|jll| . In 
tau decays leading to v e or D^, we use the distributions 
for those branches Jl^,|ll[. There are many tau decay 
branches that yield charged pions or kaons. However, 
these mesons will be stopped before they decay, so that 
the produced neutrinos are at low energies and hence ig- 
nored. For all particles, we impose a low-energy cutoff at 
100 GeV, though the results are insensitive to the cutoff 
as long as it is well below the transparency energy. 

In Fig. U (cf . Fig. 2 of Ref. ) , we plot the numerical 
results for the emergent and v r + fluxes using sev- 
eral initial energies. A very small fraction (about 1%) of 
v r either downscatter to very low energies by repeated 
NC interactions, or the tau lepton emerges from Earth; 
these are neglected in Fig. |^. The results are in general 
agreement with the analytic estimates, though the 
flux is about twice as large as the analytic estimate. The 
principal reason for this is the fact that the secondary 



3 




30 40 50 60 
■^(degrees) 



90 



FIG. 2. The surviving fluxes associated with the initial v T 
and ffj, fluxes versus nadir angle, for a variety of energies. The 
results were calculated with our Monte Carlo code, and the 
slight jaggedness in the curves is due to statistical fluctua- 
tions. The initial flux is depleted (v e would be the same), 
and the flux associated with the initial v T is increased due to 
the secondary Z/ M (and the equal i/ e component, not shown). 
The initial energies are 10 5 , 10 6 , 10 7 , 10 8 , and 10 9 GeV, go- 
ing from top to bottom for the i/ M fluxes and bottom to top 
for the v T + z/ M fluxes. The emergent v T flux is unity for all 
energies. We use the density profile given in the Preliminary 
Reference Earth Model described in, e.g., Ref. M. Note the 
mantle-core transition at ip ~ 33°, corresponding to a radial 
distance of 3480 km. 



are about 1/2 as energetic as the regenerated v T , and 
so have a longer interaction length than assumed in the 
analytic estimate. They are created in the last few inter- 
action lengths of the v T . We find that the v T and are 
distributed around the transparency energy, as expected. 
The surviving have the initial energy at least un- 
til their flux is greatly suppressed. A very small fraction 
will be moderated down to the transparency energy by 
repeated NC interactions; this is why the flux in Fig. || 
never completely vanishes. 

Conclusions. — Earth absorbs high-energy v e and v^. 
In fact, high-energy v r are absorbed too, but the v T flux 
is regenerated by prompt tau decays, thus moderating 
the v T energies down to near the transparency energy. 
As a function of the nadir angle, the v e and fluxes 
are exponentially absorbed, while the v T flux remains un- 
changed by passage through Earth. Halzen and Saltzberg 
have proposed that these v T can be detected by CC in- 
teractions after which the tau decays to a muon below 
the detector, and that they can be separated from 
CC interactions by their characteristic nadir-angle de- 
pendence 0. A difficulty with this technique is the 18% 
branching fraction for the tau decay into a muon. 

We have pointed out a new effect, that the v T — » 
t — > v T regeneration process creates a secondary flux. 
Though their flux is at most 0.2 of the v T flux, they 
are as detectable by the production of muons as the en- 



tire v T flux. Similarly, a secondary v e flux is created 
that doubles the detectability of v T by the production of 
electrons and their associated showers. The detectabil- 
ity of v T by neutral-current channels will be about 40% 
larger. For astrophysical antineutrinos, the secondary 
neutrinos created have the same effect, again doubling 
the detectability of the v T flux. Taking these secondary 
fluxes into account in the energy and nadir-angle distri- 
butions (e.g., in studies like Refs. ]lO|jL4| )) will be essential 
for understanding the spectra and flavor composition of 
astrophysical neutrinos. 
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